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The mass gap and vacuum energy of 
the Gross-Neveu model via the 2PPI expansion 
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We introduce the 2PPI (2-point-particle-irreducible) expansion, which sums bubble graphs to all 
orders. We prove the renormalizibility of this summation. We use it on the Gross-Neveu model to 
calculate the mass gap and vacuum energy. After an optimization of the expansion, the final results 
are qualitatively good. 
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I. INTRODUCTION 



The Gross-Neveu (GN) model is plagued by in- 
frared renormalons. The origin of this problem lies in 
the fact that we perturb around an instable (zero) va- 
cuum. A remedy would be the mass generation of the 
particles, connected to a non-perturbative, lower value 
of the vacuum energy. Such a dynamical mass must be 
of a non-perturbative nature, since the GN Lagrangian 
possesses a discrete chiral symmetry. A dynamical mass 
is closely related to a nonzero vacuum expectation value 
(VEV) for a local composite operator (i.e. ipip). This 
condensate introduces a mass scale into the model. We 
consider GN because the exact mass gap and vacuum 
energy || are known. This allows a test for the reliability 
of approximative frameworks before attention is paid to 
dynamical mass generation in more complex theories like 
SV(N) Yang-Mills f§. The last few years, several meth- 
ods have been proposed to solve this problem and get 
non-perturbative information out of the model [?] . 
In this paper, we address another approach, the so- 
called 2PPI expansion. Its first appearance and use for 
analytical finite temperature research can be found in 
H [jj, 11, In Secjn[ we give a new derivation of 
the expansion. Sec. Ill is devoted to the renormalization 
of the 2PPI technique. Preliminary numer ical results, 
using the MS scheme, are presented in Secfy} We re- 
cover the N — ► oo approximation, but we encounter the 
problem that the coupling is infinite. In Secjv] we opti- 
mize the 2PPI technique. We rewrite the expansion in 
terms of a scheme and scale independent mass parameter 
M . The freedom in coupling constant renormalization is 
reduced to a single parameter bo by a reorganization of 
the series. We discuss how to fix bo- Numerical results 
can be found in Sec. VI. We also give some evidence to 
motivate why results are acceptable. We end with con- 
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elusions in Sec. VII 



II. THE 2PPI EXPANSION 

We start from the (unrenormalized) GN Lagrangian in 
two-dimensional Euclidean space time. 



£ = ^-i 9 2 (^) 2 



(1) 



This Lagrangian has a global U(N) invariance and a dis- 
crete chiral symmetry ip — > j^i/i which imposes (ipi/)) = 
perturbatively. This model is asymptotically free and has 
spontaneous chiral symmetry breaking. As such, it is a 
toy model which mimics QCD in some ways. 
First of all, we focus on the topology of vacuum diagrams. 
We can divide them in 2 disjoint classes: 

• Those diagrams falling apart in 2 separate pieces 
when 2 lines meeting at the same point x are cut. 
We call those 2-point-particle-reducible or 2PPR. x 
is named the 2PPR insertion point. FIG.l depicts 
the most simple 2PPR vacuum bubble. 




FIG. 1: 
point. 



A 2PPR vacuum bubble, x is the 2PPR insertion 



• The other type is the complement of the 2PPR 
class, we baptize such diagrams 2-point-particle- 
irreducible (2PPI) diagrams. FIG. 2 shows a 2PPI 
bubble. 

We could now remove all 2PPR bubbles from the dia- 
grammatic sum building up the vacuum energy by sum- 
ming them in an effective mass. To proceed, we must use 
a little trick. Let's define 



(ijjijj) = (tpti/ji) 



(2) 




FIG. 2: A 2PPI vacuum bubble. 






FIG. 3: Generic vacuum bubble. 



where the index i = 1 . . . 2N goes over space as well as 
internal values. Obviously, we have 



/- A 



(3) 



We now calculate where E is the vacuum energy. 
The g 2 derivative can hit a 2PPR vertex or a 2PPI 
vertex, (see FIG. 3 and FIG. 4) In the first case, we have 
diagrammatically the contribution 

(Mm - tiiSjk) A fei = -\ (\ - A 2 (4) 

In the second case, we can unambiguously subdivide the 
vacuum diagram in one maximal 2PPI part, which con- 
tains the vertex hit by ^j, and one or several 2PPR 
parts which can be deleted and replaced by an effective 
mass to. A simple diagrammatical argument gives 



mSij = -g (SijSki - SuSjk) Am 



Summarizing, we have 
dE 1 



-A 2 I 1 — 

2 \ 2N 



as 



2PPI /_ 2 \ 



as 



(5) 



(6) 



(7) 



The g 2 dependence in E 2 ppi comes from the 2PPI ver- 
tices. To integrate (Q), we use the Anzatz 



E (<? 2 ) = E 2PPI (to, g 2 ) + eg 2 A 2 



(8) 



with c a constant to be determined. Using (^|), we find 
from (0) 



d£ dE 2 PPI / 2\ 

dj 2 = ^g 2 -^ 9 ^ 
,dA A 1 
2iV 



as 



2PP7 



- A,1 -2lV 
dA 



^— (l-^))+^ 2 + 2c ff 2 A^ 



(9) 



FIG. 4: Diagrammatic depiction of 
the bubble of FIG.4. 



d 



(fat dot) applied on 



A simple diagrammatical argument gives 
dE 2 ppi / 2 \ . 



(10) 



This is a (local) gap equation, summing the bubble 
graphs into to. Using ( JlOj ) and comparing (^|) with (0), 



we find c 



so that we finally have that 



\ g2 V~ 2^) A 2 +S 2 p P/ (to, 5 2 ) (11) 







12) 



It is easy to show that the following equivalence hold 

8E 2 ppi _ ^ dE 
dm dm 

One shouldn't confuse (|l^) with the usual procedure of 
minimizing an effective potential V(tp) with respect to 
the field variable ip. First of all, to is not a field variable. 
Secondly, the expression for E in terms of the 2PPI ex- 
pansion is only correct if the gap equation is fulfilled. 



III. RENORMALIZATION OF THE 2PPI 
EXPANSION 

Up to now, we haven't paid any attention to diver- 
gences. We will now show that an equation such as ( fTl| ) 
is valid for the vacuum energy E with fully renormalizcd 
and finite quantities. Since in the original Lagrangian 
there is no mass counterterm, one could naively expect 
problems with the non-perturbative mass to, which gen- 
erates mass renormalization in E 2 ppi. Another possible 
problem is vacuum energy renormalization. Perturba- 
tively, the vacuum energy is zero and hence no vacuum 
energy renormalization is needed. Non-perturbatively, 
we expect logarithmic divergences proportional to to 2 
for E 2 ppj. As we will show, both these problems are 
solved with coupling constant renormalization. 

The trick is to separate the contribution of the coupling 

constant renormalization counterterm —^SZ^ 
into 2PPR and 2PPI parts, corresponding with the 
topology of the original divergent subgraphs. Let i and j 
be the indices carried by the lines meeting at the 2PPR 
vertex, then we have 



SZ4 (SijSki - SuSkj) 



Z 7 2PPI 



r 7 2PPR 



(13) 



XX 



FIG. 5: A diagrammatical identity. 
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FIG. 7: Divergent subgraphs of ^ containing the 2PPR 



vertex s. 



and 



FIG. 6: Divergent subgraph containing the 2PPR vertex x. 
Fat lines denote full propagators. 



Note that crossing will change a 2PPR part into a 2PPI 
part. 

Because of the diagrammatical identity shown in 
FIG. 5, with m — x — n a Vv0n insertion, we have a 
relation between the 2PPR part of coupling constant 
and mass renormalization. 



(18) 



S%$8 = fa*n 



SinSmj) SZ2 



,kl 



(14) 



This identity can be used to show that the divergent effec- 
tive mass m, given by @, gets replaced by a finite renor- 
malized mass ran = Z^rn = —g 2 (l — j^) An, where 
A R = Z2A is the finite, renormalized expectation value 
of the composite operator ipip. Indeed, let us consider a 
generic 2PPR subgraph or bubble graph with a 2PPR 
vertex x. The divergent subgraphs of this bubble graph, 
which do not contain x, can be made finite by the usual 
counterterms for wavefunction and coupling constant re- 
normalization. The resulting effective mass will be given 
by @, but now with A = (ipX) evaluated with the full 
Lagrangian, i.e. including counterterms. We still have 
to consider the subgraphs of the bubble graph which do 
contain the 2PPR vertex x. They can be made finite by 
coupling constant renormalization, but because the sub- 
graph is 2PPR at x, only the 2PPR part of the counter- 
term has to be inserted and we get the contribution (see 
FIG.6) 

—g 2 SZl^ p ^Aki — —g 2 (5ijSZ2;mm,kk — SZ2;ij,kk) 7777 
' J ' ziv 

_ (15) 

where use was made of (y|) and (jlj) . Since for a diagonal 

mass matrix, we can define 8Z2 by 



5Z 2 5ki = SZ 2 



mm.kl 



we have 



2N6Z 2 = SZ- 



2:mm.kk 



(16) 



(17) 



After substition of (E7|) and into (fil|, we find that 
the contribution of the 2PPR counterterm insertion gives 



1 

2N 



~ S *j9 2 ( 1 - — ) SZ * A 



and hence a mass renormalization 



Sm = -g<[l- — )5Z 2 A 



(19) 



(20) 



so that we obtain a finite, effective renormalized mass 

1 



m R = Z 2 m = -g 2 ( 1 - — ) A R 



(21) 



with A^ = Z2A = Z2 \tptp) the finite, renormalized 
VEV of the composite operator ^■0. 

To obtain a finite, renormalized expression for the 
vacuum energy as a function of Ar or Thr, we have to 
use the same trick as in the unrenormalized case and 
consider the renormalization of 4§ . Let us first consider 

as 

the case when the vertex x hit by is a 2PPR vertex 
and restrict ourselves to divergent subgraphs which 
contain x (the ones not containing x pose no problem 
and simply replace the original A evaluated without 
counterterms by A with counterterms included). The 
divergent subgraph can just end at x from the left 
or the right (FIG. 7a and 7b) or the 2PPR vertex x 
can be embedded in it (FIG. 7c). Graphs 7a en 7b 
can be made finite by the 2PPR part of the coupling 
constant counterterm and making use of (|l|), their 
renormalization contributes 

(7a) + (76) = 2 (^-^j A ki^Z 2 -ijM {8ij6 pq - S iq S PJ ) A pq 
= ~ (l - ^) SZ2A 2 (22) 
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FIG. 9: Coupling constant renormalization graph with 2 
2PPR vertices. 



FIG. 8: Selfenergy subgraphs needing mass renormalization. 



where we have used (||), ( |l4| ) and (|l7|). Graph 7 c can 
be made finite with that part of coupling constant re- 
normalization that factorizes at the 2PPR vertex x. Its 
renormalization therefore contributes 



(7c) 



A k l (5Z2;ii^kldZ2 ; jj, pq 

1 _ -L) 5 Z 2 A 2 

OAT 2 



I'll 



(23) 



where we made use of ([17] ) an d ([18| ) . Adding the counter- 
term contributions ( 22h and (|23|) to the original unrenor- 
we obtain -| (l - ^) (Z 2 A) 2 = 
When 



malized expression 
fl 



l 

-i (l - ^) A%, which is finite. When ^ hits a 2PPI 
vertex, we can unambiguously subdivide the vacuum di- 
agrams in a maximal IP PI part, which contains the ver- 
tex hit by -ggx, and one or more 2PPR bubble insertions 
which, after renormalization, can be replaced by the ef- 
fective renormalized mass rnp. We therefore have 



dE 
dg 2 



1 



1 

27V 



dE- 



ippi i 



dg 2 



(m R ,g 2 ) (24) 



We still have to show that the usual counterterms make 
dE Qg2 FI ijnR^g 2 ) finite. The non-perturbative mass Wir, 
running in the propagatorlines, will now generate selfen- 
ergies which require mass renormalization, which is not 
present in the original Lagrangian. Again coupling con- 
stant renormalization will solve the problem. Let us con- 
sider a generic selfenergy subgraph which needs mass re- 
normalization. Since the divergence is linear in tor, we 
can restrict ourselves to 2PPI diagrams with only one 
2PPR bubble insertion (FIG. 8). The divergent part of 
this subgraph, that one wants to renormalize, can end 
at the 2PPR vertex (FIG. 8a) or can continue through- 
out the 2PPR bubble (FIG. 86). In the first case, one 
needs the 2PPR part of coupling constant renormaliza- 
tion which contains only one 2PPR vertex (because the 
divergent part considered belongs to the 2PPI part of 
the diagram). We obviously have 



5Z 



2PPR.1 



J A\ij,kl ~ {fiijfimn ^inS m j)SZ 2 . mnk i (25) 

so that the countertcrm contribution is 



fa \ 2 A x^PPR.l 

(8a) = -g A k i8Z 4 . klij - 



7 2PPI i 



(26) 



where use was made of ( pj| ) and (|25|). 
In the second case, the divergence factorizes into a 
2PPR coupling constant renormalization part (the bub- 
ble graph part) and a 2PPI mass renormalization part, 
so that the counterterm contribution is 



(86) = -g 2 A kl SZlZ% 6ZiZ,H 
1 



-g A l 



2N 



6Z 2 5Z 2PPI 8 l 



(27) 



Adding both contributions, the relevant parts of the cou- 
pling constant counterterms give 



3a) + (86) = - fl -%All- 

— xry2PPIx 

= m R dZ 2 dij 



1 

2N 



SZ. 



2PPI } 



(28) 



which is exactly what we need for mass renormalization 
in E 2 ppi- 

In an analoguous way, we can consider the loga- 
rithmic overall divergences of the vacuum diagrams 
which are quadratic in tur. We now consider 2PPI 
vacuum diagrams with two bubble insertions. One 
type of coupling constant renormalization subgraphs 
end at both 2PPR vertices (FIG.9). They can be 
renormalized by the corresponding 2PPR part of the 
coupling constant renormalization counterterm. 



&Z 2 ^ k { 2 — g 2 (SijS mn — 6i n 8 m j) (5Cnm rs {^rs^kl ~ S r i5 ks ) 

(29) 

where 6£ m n,rs is the overal divergent part of 
(4> m 4>nip r 4's)- Adding the contributions from cou- 
pling constant renormalization graphs which also go 
through the bubble parts, we find 



dE. 



2PPI 



1 X ry2PPR,2 A.R a.R 1 2 t/2PP/ 

-oZ A .- A ui A ij A kl = -m R dC 



4;ij,kl 



R°( 



with 



<K 



2PPI 



2PPI 
mm.nn 



(30) 



(31) 



and use was made of (29). 

Again coupling constant renormalization provides us 
with the necessary additive renormalization of the 2PPI 
vacuum energy. Furthermore, completely analogous ar- 
guments can be used to show that the unrenormalized 
gap equation (Hi]) gets renormalized to 



dE. 



2PPI 



dmR 



(mR,g) 



(32) 
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FIG. 10: E 2 ppi. 



It is clear that the 2PPI coupling constant and wave 
function renormalization subgraphs can be renormalized 
with the original counterterms. We therefore conclude 



that 



dE 2 



dg< 



— ( m R,g) is finite and hence (|24| ) is finite and 



can be integrated. Making use of the gap equation (32), 
we find 

E {g 2 ) = \g 2 (l - A 2 , + E 2PPI (m R , g 2 ) (33) 

Of course, we also have the equivalence ( |L2] ) in the 
renormalized case. 



For the rest of the paper, it is implicitly under- 
stood we're working with renormalized quantities, so 
that we can drop the i?-subscripts. 



IV. PRELIMINARY RESULTS FOR THE MASS 
GAP AND VACUUM ENERGY 

FIG. 10 shows the first terms in the loop expansion for 
E2PPI- Restricting ourselves to the 1 loop vacuum bub- 
ble, we have in dimensional rcgularization with d = 2 — e, 



E- 



2PPI 



-2Nm 2 - 



1 



-/i 



d d p 



1 



2-s'~ J (2Tr) d p 2 +m 2 

Using the MS scheme, we arrive at 

1 1 ( 1 \ ,? N ? A m 2 
E = -g 2 1 ) A + — to 2 In -5- - 1 

2 J V 2N J 4?r V M 



The gap equation ( |32[ ) gives 

Nm , m 2 — 
—— In — ~- = A 
2tt JI 2 

Consequently, the vacuum energy is expressed by 
E = m 2 



47T 



At 1 loop order, we have 

tm = 



A)ln# 



(34) 



(35) 



(36) 



(37) 



(38) 



where /3q is the leading order coefficient of the /^-function 



= Pit) = -2 (Pog 4 + a x f + (3 2 f + 



(39) 



The values of the coefficients can be found in 11 3L fL4L 1 



Pi 



N -1 

N 



1 



$2 = ~ 



An 2 
(N-l) 



IN- 



16tt 3 



(40) 
(41) 
(42) 



To get a numerical value for the mass gap Q > we have to 
choose the subtraction scale ~p. The choice immediately 
coming to mind is setting Ji = m, which eliminates the 
potentially large logarithm present in (|35|). Doing so, we 
find, next to the perturbative solution m = 0, 



while 



m = m = Ajjg 



The exact mass gap is given by 
mexact = (4e)^~ 



\ x 2N-2 J 



A MS 



while the exact vacuum energy is 
1 



E r 



. cot 



2{N - 1) 



(43) 



(44) 



(45) 



(46) 



We expect that the error on E consists of the error on the 
mass squared and the error on the function multiplying 
that mass squared. Therefore we will consider the quan- 
tity \J — E to test the reliability of our results. We define 
the deviations in terms of percentage P and Q, i.e. 



P 



Q 



100- 



100- 



"eff TYlexact 
^exact 

l ~^E 



(47) 
(48) 



Looking at TABLE Q, we notice that our results [S| are 
quite acceptable. We notice there is convergence (P — *• 
and Q — * 0) to the exact result in case of N — > 00. 
In fact, we recovered the N — > 00 approximation. For 
comparison, we also displayed the next to leading results 
p2|, given by expanding (|45l) and (f46|) in powers of 1/JV. 



mi/jv 



= 1 



1 -7e +ln4 



N 

El/N - I "47 



2iV 
7s — In 4 

4-7T 



°(^)) A M5( 49 ) 



6 



TABLE I: 1 loop results for mass gap and vacuum energy 



where 7(<? 2 ) governs the scaling behaviour of A 



N 


p 




J 1 /N 


Ql/N 


2 


-46.3% 




-21.9% 




3 


-32.5% 


-6 7% 


-12.2% 


O.O /U 


4 


94 9% 


-o.u /o 


7 n% 

1 .u /o 


J- .o /o 


5 


-19.1% 


-7.2% 


-4.5% 


0.4% 


6 


-15.8% 


-6.2% 


-3.1% 


0.1% 


7 


-13.5% 


-5.5% 


-2.3% 


0.007% 


8 


-11.7% 


-4.8% 


-1.8% 


-0.03% 


9 


-10.4% 


-4.3% 


-1.4% 


-0.04% 


10 


-9.3% 


-3.9% 


-1.1% 


-0.04% 


20 


-4.6% 


-2.0% 


-0.3% 


-0.02% 



dA , 2\"a" 



with 



7(5 2 ) = 7o5 2 + 7i.9 4 + 7 2 .9 6 H 

The coefficients are given by Q [TJ, [H| 



where 7b ss 0.577216 is the Euler-Mascheroni constant. 



However, the choice \x = 



Ajrq cannot satisfy us, 



70 = 

71 = 

7 2 = 



N — - 

IV 2 



7T 

AT- 



4tt 2 



47T 3 



(53) 
(54) 

(55) 
(56) 
(57) 



since we are expanding in g (m) = oo. We may have After some calculation, we find 



qualitatively good results, but for a field theory where 
the exact results are unknown, g 2 = oo gives by no means 
an indication about how trustworthy our approximations 
are. It is clear we must find a better method to achieve 
results with the 2PPI expansion. 



V. OPTIMIZATION AND 2 LOOP 
CORRECTIONS 

A. Renormalization group equation for E 

A standard approach to get better results is the usage 
of the renormalization group equation (RGE). In our 
approach, we first solved the gap equation and then 
set Jl = m. Normally, when minimizing the effective 
potential V, one first sets Jl = m, and afterwards 
the RGE is used to sum leading logarithms, while all 
quantities are running according to their renormalization 
group equations at scale m. We already mentioned E 
cannot be treated on equal footing with an effective 
potential due to the demand that J= = must hold. 
We first point out why this also disturbs a standard 
RGE improvement of E. 

Since E is the vacuum energy, it is a physical quantity 
and therefore, it shouldn't depend on the subtraction 
scale Jl. This is expressed in a formal way by means of 
the RGE 



dE 
Jl— = 
op 



(51) 



In a perturbative series expansion, this means the dif- 
ferential equation (^l[) must be fulfilled order by order, 
when all quantities obey their running w.r.t. Jl. Out of 
( |l6| ) , we extract the running of m, namely 



dm 
dji 



Pit) 



+ j(g 2 ))m = K(g 2 )m (52) 



_dE 

dji 



1 



— 9 7^_ 2 n d 



4*1-3* 



(58) 



It seems that E doesn't obey its RGE. Perturbatively, it 
of course fullfills the RGE up to O (g 2 ) since m = to all 
orders in perturbation theory. We must not be tempted 
to interpret this failure as the need to introduce some 
non-perturbative running coupling constant, as can be 
found in literature sometimes. The nature of the apparent 
problem lies in the fact that we forgot about the gap 
equation ^= = 0, because only then our 2PPI expression 
for E is meaningful. (36) gives that 



vri 1 
In — 7t oc 
M 9 



(59) 



when we would 
we would need 
,h leading log 



li- 



lt is easy to check that (59) means that all leading 
log terms in the expansion of E are of the order unity. 
Consequently, we cannot simply show order by order 
that /If = 0. 

The problem extends to higher orders 
calculate £ up to a certain order n 
knowledge of all leading, subleading,.. 
terms. 

The above discussion reveals a possible strategy : we 
could do a (leading) log expansion for E2PP1, with a 
source J coupled to tptp. Then we could use the RGE 
for E to sum all (leading) logs in E2ppi- We leave this 
idea, because the RGE for E itself is non-linear when 
J ^ p3[ . This is accompanied with its own problems. 
A thorough discussion of this subject can be consulted 
in§. 

We conclude that we cannot use the RGE for E to 
optimize that what we did hitherto. The crucial point 
is that the gap equation must hold for consistency. We 
can only set JL = m in after deriving E w.r.t. m and 
solving this gap equation, not before. 



B. Optimization 



Next, we rewrite to in terms of M by inverting ( |60| ) 



We have seen that the MS scheme is not optimal 
for the IP PI expansion used on GN. We could have 
renormalized the coupling constant in another way and 
hope that this gives better results. It is easily verified 
that going to a scheme with coupling g 2 , determined 
at lowest order by g 2 — g 2 (l + bog 2 ), gives the same 
results as in ([43] ) and (|44|), but now with g 2 = by 1 . This 
means results are as good as before, but for a sufficiently 
large bo, g 2 is small. Again, we put ~fi — m to cancel 
logarithms. 

Till now, we kept to as the mass parameter, however 
we should go to another scheme for this quantity too. 
The results are then no longer independent of the 
renormalization prescriptions, i.e. if to = to (l + aog 2 ) 
at lowest order, then ao enters the final results, and ao 
is completely free to choose. We tackle the problem of 
freedom of renormalization of the coupling constant and 
mass parameter in 4 consecutive steps. 



Step 1 

First of all, we remove the freedom how the mass 
parameter is renormalized. We can replace to by an 
unique pl| M such that M is renormalization scale 
and scheme independent (RSSI) @. Out of ([52]), we 
immediately deduce that 



M = f(g 2 )m 



where f(g 2 ) is the solution of 



-Of - f -2s-p 



(60) 



(61) 



When we change our MRS, we have relations of the form 
f = g 2 (l + bog 2 + b l9 A + ---) (62) 
to = to (l + mog 2 + TOig 4 + • • • ) (63) 

7(f) = f(g 2 ) (i + fog 2 + fig* + ■■■) (64) 

Whenever a quantity is barred, it's understood we're 
considering MS, otherwise we're considering an arbi- 
trary MRS [ p5[ . Using the foregoing relations, it is easy 
to show the scheme independence of M. 

The explicit solution, up to the order we will need 
it, is given by 



f(f) = it) 



2.1„ 



Pilo 
" Pi 



2i 

Po 



—4 

w 



fhjo 

" /3n 2 ^ A) 



— \ 2 

7i 



W (1 



cig 2 



C23 4 ) 



(66) 



where 



ci 



§J/J0 _ 
PI _ 

PI 

Pill 
PI " 



7i 
Po 



(67) 



— \ 2 



1± 

Po 



Po 



72 

Po 



(68) 



Step 2 

Transformation (^6|) allows to rewrite E in terms of M. 
Since the next contribution to ( |35| ) is proportional to 
fm 2 (see FIG. 10), we can rewrite E up to order ~g 2 



when (Bq) is applied. Explicitly, 



E 



M 2 (g 2 f-W t 
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(1 + 2 Cl 5 2 ) 



x In 
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^2 + 2ci + (2c 2 + c\) f 



(69) 



It is important to notice that the demand J= = 
is translated into = 0, because M and to differ 
only by an overall factor / which depends solely on ~g 2 (jl) . 

Step 3 

(|69|) is still written in terms of g 2 . Using (|6^), we ex- 
change ~g 2 for g 2 , where the 6, parametrize the coupling 
constant renormalization. We find 
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M 2 (g 
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(70) 



with 



8 



ei 



/r 
°0 



- 61 + c\ + 2c 2 
2(l-w) 



+ 




— (b U + 2ciV) 

47T 

2 - ^ 
Po 
7-2 



(73) 
(74) 

(75) 



Step 4 

Consider (70). We notice that the degrees of freedom, 
concerning the scheme, are settled in the bi. When we 
rewrite the expansion in terms of g 2 loop instead of g 2 , all 
scheme dependence is reduced to one parameter, namely 
b . This was also recognized in J?j]. It is in a way more 
"natural" to rewrite a perturbative series in terms of 
9iloop> because g 2 itself is changed whenever we include 
the next loop order, while g\i oov of course remains the 
same. 

The necessary formulas are given by 
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(77) 



A is the scale parameter of the corresponding MRS. In 
|l7|, it was shown that 



A = A 



MS 



o_ 

e 2 ' 3 o 



For p 2 , we have @ 



P2 = (&g-&i)A>+j9i&o+/3s 



(78) 



(79) 



Since ( |70| ) is correct up to order g 2 M 2 , we can expand up 
to order x~ x M 2 . Using (||), @ and @, the vacuum 
energy becomes 



E = AT I - 

. a: 



. 9 10 
1\ W 



iJ-lX + So + — ) (M)) 

a; ' 



FIG. 11: Diagrams needed to calculate m e // in function of 
m. 



with 
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C. 2 loop corrections 

The next order corrections are 2 loop for the mass (the 
setting sun diagram of FIG. 11) and 3 loop for the va- 
cuum energy (the basket ball diagram of FIG. 2). We 
will restrict ourselves to 2 loop corrections. The diagram 
displayed in FIG. 11 gives a mass renormalization. The 
double line is the full propagator Sf u u(p). We first em- 
ploy the MS scheme again for the calculation. 
Let V be the value of the (amputated) setting sun dia- 
gram. Since 



S(p) 



1 



ijj + m 



(87) 
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we have 



results in 



S2loop{p) = 



1 



ijj + m — V 



The effective mass rn e ff is the pole of S2i 00 p{p)- From 
pi, we obtain 



where 
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Working up to order ]p , we find for the inverse propaga- 
tor 
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(97) 



The quantities ci, c%, U, W and L are the same as defined 
before. Again, only bo is left over as scheme parameter. 



VI. SECOND NUMERICAL RESULTS FOR THE 
MASS GAP AND VACUUM ENERGY 

We first discuss how we can fix the parameter bo in 
a reasonable, self- consistent way. A frequently used 
method is the principle of minimal sensitivity {PMS) |19]] . 
This is based on the concept that physical quantities 
should not depend on the renormalization prescriptions. 
In our case, the vacuum energy E as well as the mass gap 
m e ff are physical, so we could apply PMS. However, 
PMS doesn't always work out. Sometimes there is 
no minimum, then an alternative is picking that bo 
for which the derivative of the considered quantity is 
minimal (— > as near as possible to a minimum). Also 
fastest apparent convergence criteria (FACC) can be 
practiced. 
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Solving for the pole gives 



l eff 



N- 



16tt 2 



2 In 



2 rrf_ 

— 2 

A* 



7r 2 m 2 
— - F x + 1 - 2 In + F 2 
6 (j, 2 



(95) 



(96) 



With [i = m = A-p-g, the above equation has no sense. 

Next, we follow the same steps as executed for E 
to reexpress m e // in terms of M and x. A little algebra 



But maybe the biggest barrier to a fruitful use of 
PMS (or FACC) arises from the same origin why E 
didn't seem to obey its RGE. Just as the scale depen- 
dence of E is not cancelled order by order, the scheme 
dependence of E won't cancel order by order, so we may 
find no optimal bo, and even if we would have such bo, 
it wouldn't be certain that the corresponding E really is 
a good approximation to E exact . The same obstacle will 
arise for the mass gap rn e ff. 

Apparently, we haven't got any further. We may 
have a way out through. M, as defined in (|60|), is RSSI, 
independent of the fact that it satisfies its gap equation 
or not. The 2PPI formalism provides us with an 
equation to calculate M approximately. This equation, 
= 0, is correct up to a certain order and M is RSSI 
up to that order by construction. Hence, we can ask 
that the (non-zero) solution M has minimal dependence 
on bo. This also gives a value for bo to calculate the 
vacuum energy, because the bo for E and M must be 
equal, again because E is only correct when the gap 
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TABLE II: Optimized first order results for mass gap and 
vacuum energy (Choice I) 



N 


P 


Q 


N 
4-jyx 


2 


? 


? 


? 


3 


-22.5% 


43.9% 


0.60 


A 
-t 




zo.y /o 


u.ou 


5 


-16.8% 


17.9% 


0.55 


6 


-14.6% 


13.8% 


0.54 


7 


-12.7% 


11.3% 


0.53 


8 


-11.2% 


9.5% 


0.53 


9 


-10.1% 


8.2% 


0.53 


10 


-9.1% 


7.2% 


0.52 


20 


-4.5% 


3.4% 


0.50 



TABLE III: Optimized first order results for mass gap and 
vacuum energy (Choice II) 



N 


P 


Q 


JV 

4-KT. 


2 


? 


? 


? 


3 


19.9% 


120.7% 


0.30 


4 


4.5% 


57.2% 


0.31 


5 


0.3% 


36.8% 


0.32 


6 


-1.2% 


27.0% 


0.33 


7 


-1.9% 


21.3% 


0.33 


8 


-2.1% 


17.5% 


0.34 


9 


-2.2% 


14.9% 


0.34 


10 


-2.2% 


12.9% 


0.34 


20 


-1.6% 


5.5% 


0.35 



equation is fulfilled. Also the mass gap m e // can be 
calculated with this bo- 
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FIG. 12: M(bo) in units of A w for N = 5 (Choice I, 1st 
order) . 
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FIG. 13: ff- in units of A w for N = 5 (Choice I, 1st order). 



A. First order results 

We start from the expression ( p0| ) , but we first restrict 
ourselves to the lowest order correction. 



E 



M 



70 

da 



(E-ix + E a ) 



(98) 



Until now, we haven't said anything about the freedom in 
scale ~p. Analogously as we fixed bo , we can ask ^= = 
due to the scale independence of M. For the sake of 
simplicity, we will however make a reasonable choice for 
~p. In order to cancel logarithms, we could set 71 = M. 

We observe that In 

V 0o J x' 

we could determine ~p such that W = 0, then the danger 
of exploding logarithms is also averted. We refer to this 
as Choice II. 



We refer to this as Choice I. 
always appears in the form W = In -=2- 



TABLE and TABLE fl| summarize the corres- 
ponding results. 

Some remarks must be made. 



1) We have determined the parameter bo by re- 
quiring that 



is minimal |2§]. In FIG. 12, M(b ) is 

plotted for the case N — 5 and Choice I. FIG. 13 shows 
^jp, again for N = 5 and Choice I. The plots for Choice 
II are completely similar. 



Notice that 



c)M 



db 



is relatively small. For both choices, 



it tended to zero for growing N, f.i. 
N = 10, Choice I. 



0.045 for 



Results for the mass gap agree very well with the 
exact values for Choice II, this is quite remarkable since 
we used a lowest order approximation. Choice I gives 
almost the same results as the N — > 00 approximation. 

For the vacuum energy, the results are somewhat 
less good than those obtained with a straightforward 
A/5 calculation. 

Nevertheless, the mass gap as well as the vacuum energy 
are converging, and we retrieve the correct N — > 00 
limit. Moreover, the relevant expansion parameter 
is relatively small, and behaves more or less as a constant. 



11 



TABLE IV: Optimized second order results for mass gap and 
vacuum energy (Choice I) 



M 



N 


P 


Q 


N 


2 


? 


? 


? 


3 


-0.2% 


54.8% 


0.16 


A 
-t 


9 fi% 
-z.u /o 


oo.o /o 


n i r 


5 


-3.3% 


23.8% 


0.16 


6 


-3.7% 


18.1% 


0.16 


7 


-3.8% 


14.5% 


0.17 


8 


-3.9% 


11.9% 


0.17 


9 


-3.9% 


10.0% 


0.17 


10 


-4.0% 


8.5% 


0.17 


20 


-3.7% 


2.8% 


0.19 



TABLE V: Optimized second order results for mass gap and 
vacuum energy (Choice II) 



N 


P 


Q 


N 

4-7T.X 


2 


? 


? 


? 


3 


-4.5% 


47.7% 


0.17 


4 


-6.5% 


27.9% 


0.17 


5 


-6.1% 


19.9% 


0.17 


6 


-5.4% 


15.6% 


0.17 


7 


-4.8% 


12.8% 


0.17 


8 


-4.3% 


10.9% 


0.17 


9 


-3.9% 


9.5% 


0.17 


10 


-3.5% 


8.4% 


0.17 


20 


-1.8% 


3.9% 


0.17 
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FIG. 14: M(b ) in units of A w for N = 5 (Choice I, 2nd 
order) . 
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FIG. 15: f££ in units of A w for N = 5 (Choice I, 2nd order). 



2) For N = 2 we didn't find an optimal bo. In the 
light of the exact results ([l|]) and (fi"6|), it isn't unex- 
pected that N — 2 causes trouble. N = 2 is a maximum 
of m exact , and close to = |, which is a root of m exact . 
There is a sharp drop between 2 and |, and somewhat 
lower than |, oscillating behaviour begins. What's 
more, N — 2 and JV" = I are both roots of E exact , while 
E eX act > between them. Again there is a sharp drop 
at | with oscillation somewhat before |. Problems 
with N = 2 persist at second order too, as will be seen 
shortly. 



B. Second order results 



In TABLE IV, we present second order results for 
Choice I, while TABLE |v| displays those for Choice II. 
Just as for the first order approximation, we plotted 

M(b ) in FIG. 14, and 
Choice I. Notice that 
N = 10 



dhl 



in FIG. 15 for the case N — 5, 
is smaller at second order. For 



0.022 . Again it reaches zero for infinite 

N. Again, we weren't able to extract a value for m e // or 
E for N = 2. 



C. Interpretation of the results 

When we compare the second with the first order 
results, a strange feature immediately catches our eyes. 
For Choice I, the mass gap results are better at second 
order, while the energy results are worse. For Choice 
II, the energy results are better, while the mass gap 
performs worse (except for N = 3). To make the 
comparison more transparent, we plotted the different 
mass gap results in FIG. 16 and energy results in FIG. 17. 
One shouldn't be alarmed that second order results are 
" worse" . We see that the difference between the Choice 
I and II results at first order are relatively large, for 
m e ff as well as for E. But at second order, the results 
are almost the same for both choices, whereas -M— is 
the same. This pleases us, because these results indicate 
that the choice of ~p is getting less relevant in the final 
results at second order. The fact that both (reasonable) 
choices for the scale 77, give results that are close to each 
other and are converging to the same N — > oo limit, 
convinces us that our method is consistent and should 
give trustable results. 

Yet, there is another way to check reliability. Wc 
already said FACC could be used as an alternative to 
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FIG. 16: Different results for m, 
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Consulting FIG. 18 and FIG. 19, we are able to under- 
stand why we should have ended up with qualitatively 
good results, since 5e as well as 5 m are close to 1, even 
for small N. We also see that both choices for JI should 
give comparable results, since 5e and 5 m fit with each 
other. 



We also fixed bo by demanding that 



nimal 27 , and we found that results were less good 
than those obtained by fixing bo by means of M, except 
for small N values ^8|. However, the convergence to 
the exact results for growing N was very slow. For 
example with Choice I, Q(5) = 19.4%, Q(10) = 19.1%, 
Q(20) = 14.1%. 

An analogous story held true for m e ff, where bo was 

was minimal. 



db 



determined by demanding that 

There, the deviation from the exact results was always 
bigger |p9l , and the convergence was again rather 
slow. For example, with Choice I, P(5) = 30.1%, 
P(10) = 25.7%, P(20) = 18.5%. All this corroborates 
our conjecture that M is indeed the best quantity to fix 
bo- 

Before we formulate our conclusions, we just like 
to mention that also in case of N = 2 there exist a mass 
gap and a non-perturbative vacuum energy. We already 
pointed out why we probably didn't find an optimal bo 
with our method. The best we can do with this special 
N value, is just choosing a (physical) renormalization 
scheme, but we must realize we can easily obtain highly 
over- or underestimated values in this case and that this 
is not a self-consistent way to obtain results. 



VII. CONCLUSION 



FIG. 17: Different results for 



PMS to fix bo- More precisely, we could use a FACC on 
both the energy E as the mass gap equation = 0. 
Explicitly, define 



5k = 



Exx- 1 - E 



(99) 



measuring the relative correction of the second order on 
the first order contribution. The closer 5e is to 1, the 
better it is, as an indication that the series expansion is 
under control. The quantity 5m is defined in a similar 

fashion. Unfortunately, no bo exists such that 



or 



are zero or minimal. However, we can substitute 

our PMS results in 5e and 5m and find out what these 
give. 



This paper, which had the purpose to investigate 
the dynamical mass generation and non-perturbative 
vacuum energy of the two-dimensional Gross-Neveu field 
theory, consisted of two main parts. In the first part, 
we proved how all bubble Feynman diagrams can be 
consistently resummed up to all orders in an effective 
mass to. We showed that this m can be calculated from 
the gap equation J= = 0, whereby E is the vacuum 
energy. E is given by the sum of the 2PPI vacuum 
bubbles, calculated with the 2PPI massive propagator 
(i.e. with mass to), plus an extra term, accounting for a 
double counting ambiguity. 

We showed that the 2PPI expansion can be renormal- 
ized with the original counterterms of the model. 

A very important fact is that the 2PPI expansion 
for E is only correct if the gap equation |pl = is 
fulfilled. In this context, we discussed the renormaliza- 
tion group equation for E, and showed why E doesn't 
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FIG. 18: Se as a function of N. 
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FIG. 19: 6m as a function of N. 



obey its RGE order by order, because the requirement 
of the gap equation turns terms of different orders into 
the same order. We stress that this does not mean E 
doesn't obey its RGE, or ask for the introduction of a 
"non-perturbative" (3- function. 

To get actual values for m e ff and E, we employed 
the MS scheme, and after the classical choice ~p = m 
to cancel logarithms, we recovered the N — > oo results. 



However, the corresponding coupling constant was 
infinite, so we couldn't say anything about validity of 
the results, without the foreknowledge of exact values. 
This, combined with the uselessness of the RGE for E to 
improve calculations, compelled us to search for a more 
sophisticated way to improve the 2PPI technique. 

In the second part, we first eliminated the freedom 
in the renormalization of the 2PPI mass parameter, 
by transforming to to a renormalization scheme and 
scale independent M. The consistency relation ^| = 
was completely equivalent to ^ = 0. Secondly, we 
parametrized the coupling constant renormalization. Af- 
ter a reorganization of the series, all scheme dependence 
was reduced to a single parameter bo, equivalent to the 
choice of a certain scale parameter A. 
We fixed this bo by means of the principe of mini- 
mal sensitivity (PMS). Originally, PMS was founded 
on the logical requirement that observable physics 
cannot depend on how one chooses to renormalize. 
Translated to our case, E and m e ff shouldn't depend 
on the arbitrary parameter bo- But we showed on 
theoretical grounds why applying PMS on neither 
m e ff nor E would be valid, because analogously as E 
(m e ff) doesn't lose its scale dependence order by or- 
der, it doesn't lose its scheme dependence order by order. 

Nevertheless, we gave an outcome to the problem 
of PMS. By construction, M is scheme and scale inde- 
pendent, so we can apply PMS on this mass parameter. 
This provides us with an optimal bo to calculate M, and 
consequently E and m e ff. For the scale ~p, we made 
2 reasonable choices. These 2 choices gave acceptable 
results at first order, yet there was quite a big difference 
between them. The second order results were compa- 
rable and qualitatively good, converging to the exact 
values for growing N. 

The relevant expansion parameter was relatively 
small. We gave extra evidence why results were good, 
by using a fastest apparent convergence argument. 

We explicitly checked that using PMS on E and 
m e ff to fix bo gave worse results, and the convergence 
was very slow. 

Summarizing, we have constructed a self consistent 
method to calculate the mass gap and non-perturbative 
vacuum energy. The 2PPI expansion, as well as the 
optimization procedure, are immediately generalizable 
to other field theories. 
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